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Discrete-time Calogero-Moser Model and Lattice KP 

Equations 

F. W. Nijhoff and G. D. Pang 



ABSTRACT. We introduce an integrable time-discretized version of the classi- 
cal Calogero-Moser model, which goes to the original model in a continuum 
limit. This discrete model is obtained from pole solutions of a semi-discretized 
version of the Kadomtsev-Petviashvili equation, leading to a finite- dimensional 
symplectic mapping. Lax pair, symplectic structure and sufficient set of in- 
variants of the discrete Calogero-Moser model are constructed for both the 
rational and elliptic cases. The classical r-matrix is the same as for the con- 
tinuum model. An exact solution of the initial value problem is given for the 
rational discrete-time Calogero-Moser model. The pole-expansion and ellip- 
tic solutions of the fully discretized Kadomtsev-Petviashvili equation are also 
discussed. 



1. Introduction 



Q The classical Calogero-Moser (CM) model, [1, 2], is a one-dimensional many- 

D body system with long-range interactions, and it is described by the equations of 

motion 

m n 

(1.1) x i = g^2p\x i -x j ) , i = l,2,...,N, 

r in which p(x) = p(x\u>i,u>2) is the Weierstrass elliptic P-function with periods lu\ 

and u>2- In the degenerate cases when one or both of the periods u>\ t 2 take infinite 
values, p(x) becomes: I) x~ 2 when u)\ — ► oo and u>2 — > ioo (This is called the 
rational case, and the eq. (1.1) in which p(x) is replaced by x~ 2 will be called the 

CTn rational CM model in the following); II) (sinx)~ 2 — ^ when ui\ = 7r/2 and 0J2 — > ioo 

(the trigonometric case); III) (sinhx)~ 2 + ^ when u)\ — ► 00 and 0J2 = iir/2( the 
hyperbolic case). This classical CM model, as well as its quantum analogue, are 
Q | integrable systems, cf. the two review papers, [3], for a comprehensive account 

1} of both cases. Integrable long-range interaction models of CM type, and their 

relativistic deformations, [4] , have gained new interest recently in connection with 
the representation theory of Lie algebras and affine quantum algebras, cf. e.g. [5, 6]. 
Furthermore, there is strong evidence that the CM model plays an important role in 
understanding the universal behaviour in quantum chaos, [7]. As other longe-range 
integrable models with similar connections, we can mention the Haldane-Shastry 
SU(N) chain, cf. [8], the Gaudin model, [9], and the Gamier system, [10]. 
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In a recent article, [11], we introduced an exact (time-)discretization of the 
classical CM model in the rational case. We know that many of the conventional 
continuous integrable models have exactly integrable discrete-time counterparts 
(see e.g. some of the other contributions at this Workshop), although it is often 
highly nontrivial to obtain them. These exact time-discretizations have a great 
potential for applications (numerical integration schemes, construction of cellular 
automata, discrete "differential" geometry, physical modelling). They provide exact 
interpolations and thereby finite algorithms for the calculation of classical orbits 
or -in the quantum regime- of quantum-mechanical propagators. As such, having 
exact discrete-time analogue of the CM model could be of use in the random matrix 
approach to quantum chaos. From a more general point of view, such a discrete- 
time long-range model might exhibit some new features, from which much can still 
be learnt about the true nature of integrability in the discrete regime. 

In this paper we will review the construction of the rational (integrable) discrete- 
time CM model and also generalize this discrete model to the elliptic case, namely 
we will construct an elliptic (integrable) discrete-time CM model, which goes ex- 
actly to (1.1) in a continuum limit. The pole-expansion and elliptic solutions of the 
lattice Kadomtsev-Petviashvili (KP) equation will also be discussed in this paper. 

2. The rational discrete-time CM model 

The method used here to obtain the discrete-time CM model is based on the 
observation of [12], cf. also [13], that the dynamics of the poles of special solutions 
of integrable nonlinear evolution equations is connected with integrable systems of 
particles on the line. The connection between the pole solutions of the (continuum) 
KP equation and the CM system was found by Krichever in [14], cf. also [15]. 
Here we will perform a similar construction for the discrete case. We will first show 
that the rational pole solutions of a semi-discretized version of the KP equation is 
connected with a time-discretized version of the rational CM model. Later in this 
paper we will generalize this result to the elliptic case and we will also discuss on 
a fully-discretized version of the KP equation, i.e., the lattice KP equation. 

The semi-discretized version of the KP equation reads 

(2.1) (p — q + it — w)j = (p — q + it — u)(u + u — it — u), 

where p and q are two (lattice) parameters, u is the (classical) field, the ~ denotes 
the discrete time-shift corresponding to a translation in the "time" direction while 
the ' denotes a shift or translation in the "spatial" direction. In (2.1) the time and 
spatial variables are discrete while the third variable £ is continuous. Eq. (2.1) 
is obtained from the fully-discretized version of the KP equation of [16] (see also 

(5.1) below) by letting one of three lattice parameters tend to zero. 

It can be easily proved that the compatibility condition of the two equations 

(2.2) </> = <f> ( + (p+u-u)<f>, 

(2.3) 4> = <j>z + {q + u-u)<j>, 

leads to (2.1), and (2.2-2.3) is called the Lax representation of (2.1). 
Direct calculation shows that 



(2.4) 



1 n m 

u = — with x = 1 h 

x PI 
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is a simple rational (pole) solution to the semi-discretized version of KP equation 
(2.1). Here n(m) is an integer, which is increased by value 1 when operated by ~ 
("). If we substitute this solution (2.4) into the i/'s in (2.2-2.3) , then we find that 

(2.5) <f> = 1 1- ^-)(p+k) n (q + k) m ex P (kO 

\ kx J 

satisfies (2.2-2.3), where A; is a spectral parameter. 

Enlightened by the above simple exercise, we now suppose that 



N 1 

« 2 - 6 > " = ■„ 

8 = 1 S % 

(2.7) <j> = (l-l^jh-](p + k) n (q + k) m exp(kO, 



where Xi and 6; are independent of £ (but they depend on the time variable), and 
we are to find the conditions these Xi and 6; should satisfy such that equation (2.2) 
is valid. Substituting (2.6-2.7) into (2.2) and equating to zero the coefficients of 
(£ — x i)~ 1 an d (£ — Xi) -1 , we obtain the following equations: 

(2.8) (p+k)b = ke + Lb, 

(2.9) (p+k)b = ke + Mb, 

where the vectors b = (&i, &2, &w) T and e = (1, 1, 1, 1) T and the matrices 



N 

(2.11) M = J2 



l X i X j i,i=i X i X 3 



N 



In the above equations, b and Xi denotes the discrete time-shift as stated before, 
and the elements of the N x N matrices e 8 j are defined by (eij)ki = 6ik&ji ■ The 
compatibility of (2.8-2.9) leads to the equation 

(2.12) (LM — ML)b + k(L — M)e = 0. 

(2.12) is a discrete inhomogeneous Lax's equation. It can be readily checked that 
the resulting equations of (2.12), i.e. 

(2.13) LM = ML 
and 



(2.14) 



(L - M)e = 0, 



4 



F. W. NIJHOFF AND G. D. PANG 



are consistent and give the same discrete equations of motion of a N— particle 
system: 



N 



t 



(2.15) + ^ + £ 



0, 



V 



i=\,2,...,N, 

where Xi denotes the discrete time-shift in the reverse direction to the one of Xi . We 

will call the model, for which eq. (2.15) are the equations of motion, the rational 
discrete-time CM model. We will show below that in a continuum limit these 
equations go to that of the original (rational) CM model. 

The Lax pair for the rational discrete-time CM model is given by L and M in 
(2.10-2.11). As a consequence of eq. (2.13) we have 

(2.16) I k = I k 

for any k = 1, 2, where I]. = Tr(L k ) . Thus the construction of invariants of the 
discrete-time flow given by (2.15) is similar to the continuous case. 

3. Integrability, integration and continuum limit 

3.1. Complete integrability. In order to establish the complete integrability 
of the rational discrete CM model (2.15), we need first to establish an appropriate 
symplectic structure for the N-particle system. We start by noting that eq. (2.15) 
can actually be obtained from the variation of a discrete action, given by 



N N 



t 

(3.1) s = ^2 c{x,x) = ^2 - ^2 i °g\ x i ^2 i °g\ x i - x j\ 

\ i*j 

in which the sum over n denotes the sum over all discrete-time iterates. The discrete 
Euler-Lagrange equations 

(3.2) |£ + ^ = , i=l,2,..,N, 

yield the equations of motion (2.15). Thus the system (2.15) defines a Lagrangian 
correspondence, (or multi-valued mapping), in the sense of Veselov, [17]. Choosing 
canonical momenta via the relations 

~ dC dC 
(3-3) pi = — , pi = 

oxi oxi 

we can introduce the standard symplectic form SI = X2j=i ^Pj A dxj on the phase 
space M with canonical coordinates {xi, pi}i = i ^ . Clearly, SI is preserved under 
the correspondence (3.3), SI = SI. 
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The above discussions allow us to make the following choice of canonical vari- 
ables (Pi, Xi) of the ith particle 



N 1 

(3.4) Pi=pi + ^2 



i=i 

leading to the standard Poisson brackets 

(3.5) {P i; xj} = 6ij, {Pi, Pj} = {x i; xj} = . 

It follows that in terms of these canonical variables, the Lax matrix L can be written 
as: 

N N 



(3.6) L = J2Pien-J2 



i = l ',J = 1 1 3 

Eq. (3.6) is the usual _L-matrix for the (continnum) rational CM model, so the 
fundamental Poisson bracket structure between L and L is the same as in the 
continuum case, i.e., 

(3.7) {L ® L} = [7*12, L ® I] - [r 2 i, I® L], 

where the r-matrix is given in [18]. 

Thus, in going from the continuous to the discrete CM model, both the L- 
operator as well as the r-matrix remain the same. What changes is the M-matrix, 
which will now depend on the canonical variables only implicitly. 

The involutivity of the invariants 

(3.8) {I k ,Ii} = {Tr(L k ),Tr(L , )} = for all k, I = 1, 2, .... , 

that follows as a consequence of the r-matrix structure (3.7), will lead to the in- 
tegrability of the discrete-time model by an argument presented elegantly in [19], 
forming the discrete counterpart of the Arnol'd-Liouville theorem, [17]. Therefore, 
we have proved the complete integrability of the discrete-time CM model. 

3.2. Exact integration. The integration of the rational discrete CM model 
(2.15) is performed roughly along the same line as for the continuum model. In the 
discrete case the initial value problem is posed by asking to find {xi(n)} for given 
initial data {xi(0)} and {xi(l) = 2^(0)}, where Xi(n) denotes the position of the 
particles Xi at the nth time-step. This is done as follows. Introducing the N x N 
matrices 

N N 

(3.9) X = ^Xieu , E = ^ e { j , 

i = l 2J=1 

and using the expressions for L and M in (2.10-2.11), it is readily shown that 

(3.10) (L-M)E = 0, 

(3.11) E(L-M) = 0, 

(3.12) XM-MX = -E , 

(3.13) XL-LX = 1-E, 
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in addition to (2.13). Now from the Lax equation (2.13) it is clear that we can put 
M = UU~ X and L = UAU -1 , where U is an invertible N x N matrix, and where 
the matrix A is constant, A = A , as a consequence of (2.13). Then, from (2.13) 
and (3.10-3.13) we obtain 

(3.14) Y = A~ 1 (Y - A _1 )A. 

where Y = U~ 1 XU . Noting the invariance of A under the discrete-time shift 
we can easily solve eq. (3.14). Further using the relation A = U(0)~ 1 L(0)U(0), 
which we calculate from the initial data, we then obtain the following result: the 
position coordinates {xi(n)} of the particles at the discrete time n, evolving under 
the discrete equation of motion (2.1b), are given by the eigenvalues of the matrix 

(3.15) Y(n) = C/(0)" 1 i(0)-" [X(0) - nL(O)- 1 ] L(0) n U(0) . 

Comparing eq. (3.15) with a similar formula for the continuous case, [3], we have 
to think of _L(0) _1 as playing the role of the resolvant, noting that we can always 
extract a constant parameter from L(0) that will become small in the continuum 
limit. Thus, eq. (3.15) shows that the discrete step Y(n) h- ► Y(n + 1) is actually 
a time-step 1 update of the continuous flow 1 . ^From (3.15) we see that at each 
discrete-time value the positions of the particles is uniquely determined up to a 
permutation of the particles. One can fix an order of the particles (thus fixing a 
branch of the correspondence) by imposing-as is usually done in the continuum 
case- that xi(n) < X2(n) < • • • < x^(n) . 

3.3. Factorization problem. We can think about the discrete-time CM model 
from the view point of a factorization problem. First, we note that the matrix M 
given in (2.11) is actually a Cauchy matrix. Thus, as a consequence of Cauchy's 
identity, [21], its determinants and its minors can be expressed as products of the 
form 

(3.16) det(M) = n»</(** -*/)(*/ - 

W k /x k - x t ) 

Furthermore, we see that from (2.13) and (3.10-3.13), we can obtain 

(3.17) {LM~ V )X — X(LM _1 ) = -M -1 , 
which tell us that 

(3.18) (LM-% = 

whose entries can be expressed as products as a consequence of the Cauchy's iden- 
tity. Then, we can think of the discrete-time flow as coming from a factorization 
in terms of a Cauchy matrix M and a "Cauchy-like" matrix K = LM -1 , whose 
entries factorize into products. Thus we can write down a dressing chain as follows: 

L = KM ^L = MK = KM ^ ... . 



1 We note here that the relevant expression given in our previous paper, i.e., the eq. (33) in 
[11], is equivalent to the above expression (3.15), which can be proved by using (3.13). A similar 
expression to eq. (3.15) can be obtained for the mastersymmetries of the discrete model, cf. [20]. 
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Furthermore, the matrix M can be indeed regarded as a dressing matrix, as it can 
be shown that multiple iterates of the form 

M(n + m)M(n + m - 1) . . . M(n + l)M(n) 

depend only on the initial and final time- variables n + m and n, but not on the 
values of the time- variable in the intermediate steps. 

3.4. Continuum limit. In a continuum limit (2.15) goes to that of the orig- 
inal (rational) CM model, which can be seen as follows. We first set 

(3.19) x i= Zi+nA, i = 1,2, ...,N, 
where A is a small (constant) parameter and 

(3.20) xi = ^ + (n + 1)A, i = l,2,...,N. 
Then, in the continuum limit, we write 

(3.21) zi = Zi + ezi + l/2e 2 ii + 

for all i = 1,2, ..,N, where £i = -j^Zi and e is the time-step parameter which, we 
suppose, is in the order of 0(A 2 ). Substituting (3.19)- (3.21) into (2.15), we get as 
the leading order term of (2.15), i.e. 

N 1 

(3.22) i; 8= - 2g V _ z)3 , i = l,2,...,N, 

j=1 \ z i z 3) 

where g = A 4 /e 2 . It is clear that (3.22) are exactly the equations of motion of 
the (continuous) rational CM model. It is interesting to note that the coupling 
constant of the continuous model arises from the discrete-time step. 

4. Elliptic generalization 

4.1. Elliptic discrete-time CM model. Now we generlize the above result 
to the elliptic case, and we are to construct an (integrable) elliptic discrete-time 
CM model, which goes exactly to (1.1) in a continuum limit. In order to do that, 
let us first recall the basic definitions and some properties of the classical functions 
of Weierstrass, see e.g. [22]. Let u>\ and u>2 be a pair of periods. The Weierstrass 
sigma-function is an entire function defined by 



(4.1) a(x) = x TT (1 - — ) exp 



— + — ) 



with uiki = kuii + lu>2- The remaining functions can be defined by the relations 

(4.2) CW = 4t ' pW = -C'W- 

Note that both the zeta-function and the sigma-function are odd functions and they 
satisfy the following fundamental relation 

4.3 C A +CG» +C" -C A + ^ + ^ = \ Y 1 \ - 

(t(\)(t{}i)(t{v)(t{\ + n + v) 
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Now by using (4.3), it can be directly checked that 
(4.4) u = C( x ) with x = £ + na + m/3 

is a simple elliptic solution to the semi-discretized version of KP equation (2.1), 
where p = ((a) and q = ((/3), corresponding to the following solution of the linear 
equation (2.2) 

a(x — k) f a(a + k)\" ( a(fi + k) \ m ^( K )^ 



(4.5) <p(x; K ) - a{x)a{K) y a{a)(T{K) ) y a{p)(T{K) ) e 

Similar to the rational case, we suppose that 

N 

(4.6) u = 

8 = 1 

N 

(4.7) <p = -*i)e <(,6)e > 

8 = 1 

where 

a(x)a(K) 

Substituting (4.6-4.7) into (2.2) and equating to zero the coefficients of $ K (^ — Xi) 
and $ K (^ — Xi), we obtain 

(4.9) pb = Lb, 

(4.10) b = Mb, 
where the vector b = (&i , b^, &jv) t and the matrices 

N ( N N 

(4.11) L ^2(( x i ^2(( x i ~ x j) e a - ^2 ~ x i) 6i i> 



8 = 1 \ ?' = 1 i=i / *,i=i 

\ i*< / <*i 

N 

(4.12) M = ^4 x i- x i) e ij- 
«' j=i 

The compatibility of (4.9-4.10), i.e., LM = ML give the equations of motion of the 
elliptic discrete-time CM model: 

(4.13) 

((xi - Xi) + ((xi - Xi) + ^ ^C(*8' - Xj) + ((xi - xj) - 2((xi - Xj)^j = 0, 

i=l,2,...,N. 

Eq. (4.13) defines an elliptic correspondence, which below is shown to be integrable 
along the same lines as the rational CM model. Apart from mappings that can be 
inferred from a lattice version of the Landau-Lifshitz equations that was derived 
some years ago in [23], this is one of the first examples of an integrable elliptic 
mapping. We remark that the elliptic case was curiously absent in the list of 
integrable mappings derived by Suris in [24]. Note that in the degenerate case 
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when u>i — ► oo and u>2 — > ioo, we recover from eq. (4.13) the rational discrete-time 
CM model (2.15). Obviously, in the intermediate degenerate limits with only one of 
the periods ui^ taking an infinite value, it goes to the trigonometric and hyperbolic 
cases. 

4.2. Complete integrability. The complete integrabilty of the elliptic discrete- 
time CM model (4.13) can be proved in a similar way to that of the rational case. 
We only note here that in the elliptic case the action is 

(4.14) 

^ N N 



\ 



',3 = 1 
<±3 



and the canonical momenta Pi in the standard Poisson brackets (3.5) now is given 

by 

dC N N N 

(4.15) Pi = -Q^ + ~ x i) = Z^to ~ x i)~ ~ x i)- 

1 3 = 1 7=1 3 = 1 

3*' 3*' 

In terms of the canonical variables, the Lax matrix L in (4.11) has again the same 
form as that of the continuum (elliptic) CM model. Thus we can use the r-matrix 
structure given in [25] to establish the involutivity of the invariants, i.e., 

(4.16) {I k ,Ii} = {Tr(L k ),Tr(L , )} = for all k, I = 1, 2, .... , 

where L is given by (4.11). 

4.3. Continuum limit. It can be shown that in the continuum limit de- 
scribed in section 3.4 the leading order term of (4.13) is exactly (1.1), where the g 



equal to -A 2 /[e 2 C'(A)]. 



5. Lattice KP equation 



In the previous sections we use pole-expansion for a semi-continuous version 
of the KP equation. In this section we will show that similar pole-expansions can 
be also done for the fully discretized KP equation, i.e., the lattice KP equation 
presented in [16], and the discrete-time CM model can be also obtained from pole- 
solutions of the lattice KP equation. 

The lattice KP equation reads 

p — r + it' — u q — r + u' — u 

(5.1) — = r , 

p — r + u' — u q — r + u' — u 

where r is the third (lattice) parameter in addition to p and q. The ' in (5.1) 
denotes a third shift or translation. Note that the three lattice directions are on 
the same footing. Eq. (5.1) is the compatibility condition of the following three 
equations, forming - as it were - an inhomogeneous linear system for (5.1), 

(5.2) <p = (p-u)(p + i), 

(5.3) (p = (q-u)(p + i), 

(5.4) ip = (r — u')ip + ip, 
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in which ip is any function independent of the particular discrete direction. Thus 
eq. (5.2-5.4) can be considered to be an inhomogeneous Lax representation of (5.1). 
Again, direct calculation shows that 

(5.5) u = C( x ) with x = £ + no. + m/3 + hj 

is a simple elliptic solution to the lattice KP equation (5.1), where p = C( a )> 
q = C(P) an d r = C(t); corresponding to the following solution of the Lax 'triple' 

a(x - k) ( a(a + k) \" ( a(j3 + k) \™ / 17(7 + k) x h 

(5.6) ip(x;K) — 1 11 11 



<t(;k)<7(k) \(r(a)(T(K) J \a(/3)a(K) J \a(j)a(K) 
ip{x;n) = [((k) + ((x - k)](p(x;k) , 

where /i is an integer similar to n and m and increased by value 1 when operated 
by '. Here and hereafter we always assume £ is a dummy variable. 
Similarly if we suppose that u and <p are given by (4.6-4.7), and 

N 



(5.7) 4> = ]>>[CM + C(e-^-K)]M€-^)e C(K)e 

8 = 1 

N 



',3 = 1 



then from (5.2) we can obtain again the elliptic discrete-time CM model (4.13) and 
its Lax representation (4.9-4.10). 

Of course, from the other two linear equations (5.3-5.4), we can get equations 
similar to (4.13), which represent the other two flows in the other two lattice direc- 
tions. These three flows are, in fact, all compatible if 

N N 

(5.8) p-q = ^2((xi - Xj) - ^((xi - Xj), 

i=i i=i 

N N 

(5.9) q-r = ^((xi - xj) - ^((xi - xj'), 

i=i i=i 

for all i = 1,2, ...,N. More details about these compatible conditions and elliptic 
solutions of the lattice KP equations will be published in [26]. 



6. Discussion 

In this paper, we have constructed a discrete-time CM model. This discrete 
model is also integrable and is the iterate of a symplectic correspondence. The 
original CM model is a limit of this discrete-time CM model, and in this limit the 
coupling constant for the long-range interaction term is encoded in the discrete-time 
step parameter. The result of the present paper fills a gap in a series of various 
deformations and generalizations of the CM model (such as the quantization, the 
deformation to the relativistic case or the generalization to different root systems), 
and it invites a number of interesting problems to be studied. First of all, one may 
ask whether one can find time-discretizations of the relativistic version of the CM 
model, [4], cf. also [27]. In fact, Suris, in [28], has found an interesting connection 
between the discrete-time Toda model and its relativistic version. Furthermore, 
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a direct connection exists between the sine-Gordon soliton solutions and the rela- 
tivistic CM model, cf. [29, 4]. This could lead to an alternative way to discretize 
the model. Secondly, for the discrete-time model, following the similarities with 
the structure of the integrable quantum mappings studied in [30], one should in- 
vestigate not only the L-part of the Lax pair, but also take the M-part under 
consideration in the r-matrix structure. This has been done for the integration of 
mappings of KdV type in [31]. Finally, these investigations should also be pursued 
on the quantum level. Although, the discrete-time model has an obvious quantum 
counterpart, in much the same way as the continuum model, the work of [30] on 
quantum mappings, as well as the work [32] on the quantum CM model, indicate 
that some important modification (e.g. with respect to the construction of exact 
quantum invariants of the discrete-time flow) might be expected. The quantum 
model in itself might shed some new light on the random matrix approach to quan- 
tum chaos, as -in principle- it provides us with a finite and exact algorithm to 
calculate multi-particle propagators in the long-range model, in the spirit of [33]. 
All these problems are being investigated and will be discussed in detail in future 
publications. 
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